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Vertical  Transport  by  Small  Seal* 
Stratospheric  Turbulence:  A  Critical  Review 


1.  INTRODUCTION 

The  meaning  of  the  word  " turbulence"  Is  ambiguous,  on  one  hand,  It  could 
refer  to  large  scale  synoptic  motions  that  take  place  on  a  global  scale,  In  this 
connection  the  phrase  "two  dimensional  turbulence"  has  been  used.  On  the  other 
hand,  It  more  frequently  refers  to  small  acale  three-dimensional  chaotic  motion 
which  causes  Intimate  mixing  on  a  small  scale.  Similarly  the  term  "eddy  dif¬ 
fusion  coefficient, "  which  Implies  a  pseudo-diffusion  effect  due  to  the  eddy  flow, 
can  be  used  In  more  than  one  manner.  On  one  hand,  It  can  Include  large  scale 
synoptic  effects  (which  are  most  often  regarded  as  advectlve  In  nature)  together 
with  small  scale  turbulence  effectB.  On  the  other  hand,  It  can  refer  exclusively 
to  small  scale  three-dimensional  turbulence  effects,  1  As  can  be  seen  by  the  title 
of  this  review,  only  the  latter  type  of  "diffusion"  will  be  considered  here. 

The  stratosphere,  by  definition,  Is  an  exceptionally  stable  part  of  the  earth's 
atmosphere,  Turbulence  In  such  a  stable  fluid  has  a  certain  peculiarity  of  struc¬ 
ture  which  must  not  be  Ignored.  It  occurs  In  relatively  thin  layers  separated  by 
what  are  usually  large  layers  of  essentially  laminar  flow.  This  Is  true  not  only 

(Received  for  publication  17  February  1901) 

1.  Tennekes,  II.  and  Lumley,  J.  L,  (1072)  A  First  Course  In  Turbulence.  The 
MIT  Press,  Cambridge,  MA.  *  “ " *  — —— — 
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In  addition,  it  is  also 


for  the  stratosphere  but  for  the  upper  ocean  as  well. 
sometimes  true  for  the  troposphere. 

The  word  "bllnl"  (pancakes)  is  sometimes  used  to  describe  thlB  layered 
structure  of  stratified  turbulence  In  the  atmosphere.  5  Such  clear  air  turbulence 
(CAT)  pancakes  are  of  the  order  of  100  m  thick  and  10  km  in  the  horizontal  dtrec 
tion,  They  are  due  to  the  shear  or  Kelvln-Helmholz  Instability.  These  layers 
are  always  assumed  to  occur  at  random  heights  and  times  with  random  thickness 
es.  They  are  not  only  intermittent  but  presumably,  also  rare  (of  the  order  of 
one  percent  of  the  fluid  volume), 

Let  Kg  denote  the  effective  diffuslvity  for  stratified  (small  scale)  transport 
over  large  regions  of  altitude.  The  subscript  B  stands  for  "bulk"  as  opposed  to 
the  local  eddy  diffuslvity  that  might  be  encountered  within  an  active  turbulent 
layer.  This  parameter  Is  the  one  upon  which  we  will  focus  our  attention.  The 
two  main  questions  which  need  to  be  answered  arei  (1)  How  Important  is  the 
role  played  by  small  scale  turbulence  in  the  overall  vertical  transport  of  tracers 
in  the  stratospherej  and  (2)  What  is  its  apptoxlmate  value? 

The  practical  importance  of  these  questions  derives,  of  course,  from  the 
stratospheric  pollution  problem.  As  is  well  known,  the  possibility  exists  that 
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Woods,  J.D.  (1008)  Wave- induced  shear  instability  in  the  .summer  thermo- 
cline,  J.  Fluid  Mech.  32(791-800, 

Woods,  J.D.  and  Wiley,  R,  L.  (1972)  Billow  turbulence  and  ocean  micro- 
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oxides  of  nitrogen  (NOx)  and  chlorofluorocarbons  can  catalyttcally  destroy  ozone 

and,  hence,  cause  global  problems  due  to  the  resulting  enhanced  transmission  of 

solar  ultraviolet  radiation  to  the  earth's  surface,  (Some  of  the  possible  disastrous 

19 

consequences  are  described  in  CLAP  Monograph  1,  1875,) 

14 

As  Reiter  has  pointed  out,  there  are  several  mechanisms  for  vertical 
transport  associated  with  the  stratosphere.  His  paper,  which  concerns  the 
stratospheric -tropospheric  exchange  processes  (not  necessarily  transport 
throughout  the  volume  of  the  stratosphere)  lists  the  following  In  order  of  Impor¬ 
tance]  (1)  Hadley  cell  motion  (38  percent),  (2)  large  scale  eddies  of  the  scale 
of  cyclones  and  anticyclones  (20  percent),  and  (3)  seasonal  tropospheric  height 
changes  (10  percent).  The  percents  refer  to  the  fraction  of  the  mass  equivalent 
to  one  hemisphere  stratosphere  transferred  through  the  tropopause  In  a  year. 

These  all  add  up  to  values  which  are  consistent  with  "residence  times"  of  mate¬ 
rials  deposited  In  the  stratosphere  such  as  atomic  bomb  debrlB.  The  "fallout 

times"  or  "resident  times,"  of  course,  measure  the  effects  of  all  processes 

14 

simultaneously,  Reiter  also  states  that,  with  regard  to  stratospheric-tropo¬ 
spheric  exchange,  turbulence  plays  an  apparently  Insignificant  role.  Whether  or 
not  such  a  conclusion  Ib  correct,  It  does  not  answer  the  question  of  whether  or 
not  small  scale  turbulence  within  the  volume  of  the  stratosphere  Is  significant 
for  vertical  transport  there. 

This  raises  an  Interesting  question]  If  fallout  times  are  known  to  some  de¬ 
gree,  why  Is  It  ImpoJtant  to  know  the  details  of  the  removal  process  and,  In 
particular,  the  relative  significance  of  the  role  of  turbulence  ?  One  answer  Is 
that  one  must  understand  the  mechanism  of  transport  If  one  Is  ever  to  estimate 
the  effects  of  large  perturbations  In  stratospheric  composition,  These  could  be 
caused  by  long-term  gradual  changes  due  to  pollution  or  by  short-term  catastrophic 
changes  due  to  large  nuclear  effects  or  rare,  natural,  large  perturbations.  Such 
large  chemical  changes  could  completely  change  the  dynamics  of  the  stratosphere 
and,  hence,  the  residence  times,  After  all,  In  spite  of  the  fact  that  the  strato¬ 
sphere  Is  defined  in  terms  of  Its  dynamic  stability,  It  Is  in  fact  Its  composition 
(that  is,  the  ozone)  which  causes  the  stability  to  exist.  The  ozone,  by  absorbing 
ultraviolet  light,  causes  the  temperature  Inversion  which  Is,  In  effect,  the  strato¬ 
sphere,  This  stability,  In  turn,  enables  large  amounts  of  ozone  to  accumulate 
without  too  much  loss.  The  ozone,  In  effect,  has  created  Its  own  container,  Its 
depletion  would  also  deplete  the  container!  Thus,  composition  affects  dynamics, 
and  a  model  which  takes  all  the  Important  mechanisms  and  these  mutual 

13.  CIAP  (1073)  The  Natural  Stratosphere  of  1074  CLAP  Monocraoh  1 

DOT  -  TS  T  -  V  5  -  3 1 .  - — - 

14,  Rotter,  E,H.  (1973)  Stratospheric-tropospheric  exchange  processes 

Rev,  Oeophys,  and  Space  Phys.  _13:45f)-474, 
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Interactions  into  account  would  bo  needed  to  estimate  the  impact  and  subsequent 
effects  following  a  large  perturbation. 

There  is  a  second  reason  for  the  necessity  of  knowing  the  value  of  Kg,  This 
relates  to  stratospheric  chemistry  In  general.  Chemicals  cannot  react  until  they 
are  mixed  Into  Intimate  contact.  Turbulence  of  the  small  scale  variety  is  the  only 
mechanism  that  con  bring  this  about,  It  Is  well  known  that  the  vertical  profile  of 
any  stratospheric  constituent  Is  highly  jagged  and  layered,  This  Indicates  that 
vertical  mixing  takes  place  over  small  vertical  scales  and  also  that  It  Is  erratic. 

One  of  the  most  Important  objectives  of  this  review  Is  to  reveal  the  high 
level  of  Ignorance  which  surrounds  the  value  of  Kg,  the  effective  turbulence  dif¬ 
fusion  parameter.  Unfortunately,  the  magnitude  of  Ignorance  Is  not  generally 

i  e  iq 

appreciated.  The  majority  of  writers  ~  seem  to  regard  the  Issue  as  relatively 
well  In  hand  and  that  the  order  of  magnitude  of  K-  (globally  averaged)  la  0. 01 
m  /s.  This  Is  to  be  compared  to  the  overall  "eddy  dlffusivlty"  based  on  fallout 

IB  •)  1 

times  and  due  to  all  processes  which  are  10  to  100  times  larger.  It  will  be 

shown,  however,  that  the  estimate  of  Kg  —  0,01  m  /•  rests  on  very  questionable 
foundations.  While  this  does  not  mean  that  this  estimate  Is  not  correct  In  principle, 
It  does  mean  that  more  experiments  will  have  to  be  performed  before  one  can 
accept  It  as  scientifically  valid. 

Our  plan  Is  primarily  to  review  the  techniques  by  which  Kg  Is  estimated  and 
to  examine  the  assumptions  (stated  and  otherwise)  which  are  Involved.  As  will 
be  seen,  thero  are  a  number  of  unanswered  questions  that  are  thuB  exposed. 

In  more  detail  we  shall  do  the  following.  First,  we  shall  consider  the  work 

1  e  4  p 

of  Lilly,  et  al,  '  They  estimated  the  value  of  the  local  eddy  diffusion 


16.  Lilly,  D.  K. ,  Waco,  D.  E. ,  and  Adelfang,  S.I.  (1974)  Stratospheric  mixing 
estimated  from  hlgh-altltude  turbulence  measurements,  J.  Aopl.  Met, 
13)488-493.  -  - 

16.  Lilly,  D. ,  Waco,  D, ,  and  Adelfang,  S.  (1978)  Stratospheric  mixing  estimated 

from  hlgh-altltude  turbulence  measurements  by  using  energy  budget  tech¬ 
niques,  The  Natural  Stratosphere  of  1974.  CLAP  Monograph  1,  DOT- 
TST -7 5-51)6-81  to  8-90, 

17.  Heck,  W.  and  Panofsky,  H.A.  (1975)  Stratospheric  mixing  estimates  from 

heat  flux  measurements.  The  Natural  Stratosphere  of  1974.  CIAP 
Monograph  1,  DOT-TST-frS-flliB-BO  to  8-62, 

18.  Mahlman,  J.  D,  and  Moxlm,  W.J,  (1978)  Tracer  simulation  using  a  global 

general  simulation  model  ...,  J,  of  the  Atm.  Set,  35)1340-1374. 

19.  Junge,  C.E.  (1963)  Air  Chemistry  and  Radioactivity.  Academic  Press, 

New  York,  250.  ~  “ 

20.  Bauer,  E,  (1975)  Dispersion  of  tracers  in  the  atmosphere,  The  Natural 

Stratosphere  of  1974.  CIAP  Monograph  1,  Final  Report,  fioT-TST-75-51, 
6-Se  to  6-55. 

21.  Luther,  F.M,  (1975)  Large-Scale  Eddy  Diffusion,  The  Natural  Stratosphere 

of  1974,  CIAP  Monograph  1,  Final  Report,  DOT-TsiT-7  5-Sl, 
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coefficient,  Kg  (H  for  heat  transfer),  from  power  density  spectra  of  wind  fluctua¬ 
tions.  In  contrast  to  Kg,  Kg  refers  to  the  diffusion  of  a  passive  scalar  such  as 

heat  that  Is  caused  by  the  active  turbulence  within  a  bllnl.  The  wind  fluctuation 

22 

measurements  were  obtained  from  the  published  reports  of  Crooks  et  al,  which 
gave  In-sttu  measurements  obtained  by  means  of  an  Instrumented  U-2  aircraft 
flying  In  the  stratosphere.  The  spectra  were  first  used  to  estimate  e,  the  viscous 
dissipation  rate.  The  latter  was  used  to  obtain  Kg.  As  mentioned,  when  they 
then  estimated  Kg,  from  Kg,  they  obtained  0,01  m2/s  as  the  order  of  magnitude. 
Next,  we  will  consider  the  work  of  Panofsky  and  Heck,  ^  They  obtained  experi¬ 
mental  evidence  based  on  the  HICAT  data  which  supported  the  technique  of 

Lilly,  et  al  for  the  estimation  of  Kr,  from  e.  The  work  of  Zimmerman  and 
23  H  22 

Loving  Is  discussed  next.  They,  too,  obtained  e  from  the  HICAT  spectra, 

but  their  method  for  estimating  K„  from  e  differed  from  the  method  of  Panofsky 

1  R  1  7  *1 

and  Lilly.  “  Their  values  for  K-,  (and  by  Implication,  using  the  reasoning  of 

1  B  IQ 

Lilly,  et  al,  Kg)  *  were,  by  comparison,  enormous. 

After  this  a  general,  critical  discussion  Is  given  of  some  of  the  basic  assump¬ 
tions  made  by  the  above  authors,  In  this  discussion,  several  unanswered  ques- 

D 

tlons  are  raised.  The  work  of  Rosenberg  and  Dewan,  Is  then  reviewed.  They 

estimated  K-r,  -  0.  3  m2/a,  a  value  In  between  those  of  Lilly  et  al,  15,  and 

Zimmerman  et  al.  The  unanswered  questions  raised  by  this  work  are  pointed 

18 

out.  Finally,  In  the  conclusion,  the  Impact  of  the  work  of  Mahlman  et  al,  on 
global  circulation  simulations  Is  pointed  out.  Two  experiments  are  given  In  the 
text  which  are  the  most  crucial  ones  to  perform  at  this  time  In  order  to  permit, 
perhaps  for  the  first  time,  a  reasonably  dependable  estimate  of  Kg. 

2.  THE  WORK  OF  LILLY,  WACO,  AND  ADELFANG 

Project  HICAT  (High  Altitude  Critical  Atmospheric  Turbulence)  Involved 
286  flights  of  U-2  aircraft  In  the  14  to  21  km  altitude  range  of  the  stratosphere. 

The  "turbulent  patches"  Included  only  2  to  6.  2  percent  of  the  total  flight  distances 
and  were  highly  correlated  with  categories  of  terrain.  As  was  mentioned  above, 
the  power  density  spectra  obtained  from  these  data  were  used  to  obtain  the  vis¬ 
cous  dissipation  rate,  e,  and  Kg  by  these  authors. 

22,  Crooks,  W.M.,  Hobllt,  F.  M.,  and  Prophet,  D.  T,,  et  al  (1687)  Protect 

HICATt  An  Investigation  of  High  Altitude  Clear  Air  Turbulence.  Tech. 
ftept.,  see  also  flP  MTIM  (1585).  - 

23.  Zimmerman,  S.  P.  and  Loving,  N.  (1975)  Turbulent  dissipation  and  dtffuslvl- 

tles  In  the  stratosphere  using  Richardson's  technique,  The  Natural 
Stratosphere  of  lt*74.  CIAP  Monograph  1,  Final  Report^  DO'f-CIAP, 
DOT-TST-75-51. 


To  obtain  t ,  they  noted  that  the  log-log  plots  of  the  spectra  had  a  nearly  -5/3 
slope,  and  hence  they  assumed  that  this  Implied  that  they  were  observing  an 
"Inertial  range"  spectrum,  As  is  welt  known,  Kolmogoroff24  predicted  this  slope 
on  the  basts  of  a  similarity  argument,  and  It  has  subsequently  been  amply  veri¬ 
fied  experimentally,  (It  Is  frequently  observed  In  geophysical  flows,)  ThuB, 
these  authors  used  the  Kolmogoroff  relation 

^(k)  =  a  e2/3  k"5/3  (1) 

(where  4(k)  Is  the  one -dimensional  velocity  fluctuation  power  density  spectrum, 
k  the  wave  number,  and  where  a  Is  the  constant  of  order  unity  derived  from 
experiment).  Solving  Eq.  (1)  for  e  gives  Its  value  In  terms  of  the  spectrum, 

More  specifically,  they  integrated  Eq,  (1)  from  the  value  of  k  corresponding 
to  wavelengths  equal  to  (810  m)'1  or  (2,  000  ft)'1  to  k  »  «e,  Solving  the  resulting 
equation  for  e  and  using  the  fact  that  the  Integral  of  ^(k)  over  k  Is  equal  to  the 
mean  square  fluctuation  velocity,  "v®,  they  obtained 


where  k^  ■  (2»/810  m)  and  Is  the  constant  of  order  unity  which  depends  on  l 
(the  velocity  component  la  designated  by  l).  In  Eq.  (2)  the  overbar  denotes  an 
average  In  the  sense  that  one  can  use  several  spectra  derived  from  a  number  of 
traverses  through  turbulence  In  order  to  estimate  the  average  of  t  . 

In  order  to  estimate  from  such  estimates  of  e  they  made  the  following 
assumption! 


P  -  B  =  e 


(3) 


where  P  means  turbulent  energy  production,  and  B  Is  the  "up  gradient  buoyancy 
flux"  or  change  of  potential  energy  caused  by  mixing.  More  exactly 


P  ■ 


B.  zL ®l«Ll 

0 


(4) 


24.  Kolmogoroff,  A.  and  Topple r,  L.  (19(11)  Turbulence.  S.  K,  Frledlandcr.  Ed.. 
Intersclence  Publisher,  NY. 


where  primed  quantities  are  the  turbulent  deviations  from  the  averages,  g  is  the 
acceleration  of  gravity,  7  the  average  potential  temperature,  V  the  average  hori¬ 
zontal  velocity,  w1  the  vertical  velocity  fluctuation,  v'  the  horizontal  velocity 
fluctuation,  9  the  potential  temperature  fluctuation,  and  Z  Is  the  vertical  coordi¬ 
nate. 

Next,  they  used  the  definition  for  K. ^  given  by 


’■H  8Z 


and  the  corresponding  definition  for  eddy  viscosity  (or  momentum  dlffuslvlty) 


v'w'  ■  *km  U 


In  order  to  obtain 


P  ■  K 


M  (llj 


Since  the  buoyancy  frequency,  Np,  Is  defined  by 


•4 -(f)  ft 


one  obtains  from  Eqs  (7)  and  (3) 


km  s2  -  kh  nb  3  * 


where  S  ■  8V/0Z,  the  vertical  shear  of  the  average  horizontal  winds.  The  defini¬ 
tion  of  the  flux  Richardson  number,  Rf  Is 


Rj.  ■  B/P 


Hence  Eq,  (9)  can  be  written 


khnb('s;- ‘)  ■■  ■ 
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Citing  the  laboratory  work  of  Thorpe,  they  set  Rj,  =  1/4  and  arrive  at 


(11) 


K 


II 


(12) 


They  inserted  values  of  e  from  Eq,  (2)  into  Eq.  (12)  to  get  the  estimates  of  Kg  in 
their  paper. 

Since  Kjj  refers  only  to  the  local  value  of  eddy  dlffusivity  within  the  thin 
actively  turbulent  layers,  they  then  proceeded  to  estimate  Kg  from 

Kg  *  F  Kg  (13) 

where  F  is  the  fraction  of  the  vertical  dimension  occupied  by  the  turbulent  layers. 
The  value  of  F  was  estimated  by  assuming  it  to  be  equal  to  the  fraction  of  the  U-2 
flight  trajectory  that  was  "turbulent."  This  is  of  the  order  of  a  few  percent,  Tak¬ 
ing  suitable  account  of  the  variation  of  F  with  terratn  category  they  concluded  that 

2 

Kg,  averaged  on  a  global  scale,  is  0,012  m  /s  to  within  (they  claimed)  one-half 

order  of  magnitude.  This  was  estimated  from  values  of  Kg  which  ranged  from 

0.4  to  1, 14  m2/s.  This  is  the  value  cited  by  Reiter14  and  Mahlman  et  al,  18  in 

connection  with  small  scale  turbulence  effects  in  the  stratosphere  and  which  seems 

to  be  the  currently  accepted  estimate. 

There  is  an  unanswered  question  In  connection  with  the  assumption  that  the 

HICAT  spectra  are  in  the  "inertial  range"  of  length  scales.  This  assumption 

seems  to  me  to  be  in  serious  contradiction  with  certain  experimental  results  pub- 

2(i  27 

llshed  in  the  literature.  Measurements  were  made  in-sltu  by  Crane,  Cadet 

28  8  ri 

and  Barat.  (See  also  Rosenberg  and  Dewan  and  Stewart.  )  Table  1  summarizes 

)ft 

Note,  however,  that  the  U-2  flew  mostly  horizontally  whereas  F  refers  to  the 
vertical  dimension, 

25.  Thorpe,  S.A.  (1972)  Turbulence  In  Stably  Stratified  Fluids:  A  Review  of 
Laboratory  Experiments,  IUCRM  Colloquium  on  Boundary  Laver 
Meteorology.  San  Diego,  cA. 

2ti.  Crane.  R.K.  (1977)  Stratospheric  Turbulence  Analysis.  Air  Force  Geophysics 
Lab.,  Final  R ep0rT.“ATOT-TR -7f- W07,  A DTff 740. 

27.  Cadet,  D.  (1975)  Vertical  wind  shear  measurements  in  the  lower  stratosphere, 

Quart.  J,R,  Met.  Soc.  101:485-493. 

28.  Barat,  J.  (1975)  Etude  experlmentale  de  la  structure  du  champ  de  turbulence 

dans  la  moyenne  stratosphere,  C.R,  Acad,  Sc,  Paris  280(Serlal  B)iG01-O93. 
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Table  1.  Measurements  of  Turbulent  Layer  Thickness 


Method 

Authors 

Turbulent  Layer 
Thickness 

Balloon-borne 

Instruments 

Barat28 

29 

Anderson 

200  m 

240  m 

Crane28 

200  m 

Rocket  trail 
derived  R, 
number 
profile 

Rosenberg  and 

Dew an8 

80  m  X2 

U-2  vertical 
traverse 

Crane28 

50  m  (light  to  moderate) 

100  m  (all) 

Crooks  et  al22 

670  m  ±  580  m 

these  measurements  of  turbulent  layer  thicknesses,  All  indtcate  that  the  order  of 
magnitude  of  the  turbulent  layer  thicknesses  is  100  m,  This  would  also  be  the  outer 
scale  presumably.  The  largest  scales  for  the  inertial  range  would  be  of  the  order 

OQ  1 

of  1/10  of  this  or  1/100,  In  other  words,  the  largest  scales  of  the  inertial 

range  are  expected  to  be  10  m  (or  even  as  small  as  1  m),  In  contrast,  the  smallest 

scales  of  the  HICAT  spectra  are  about  50  m,  the  -5/3  Blopes  are  Been  out  to  600  m 

routinely,  and,  on  occasion,  to  beyond  10  km.  A  spectrum  of  -5/3  out  to  5  km  is  not 

not  unusual.  Such  large  scales  are  not  only  Incompatible  with  the  Inertial  range 

assumption,  but  the  assumption  of  turbulent  motion  as  well.  Let  us  assume  that 

2  R  27  9  B 

the  observations  of  Crane,  Cadet,  and  Barat  are  valid  and  not  misleading  in 

any  fundamental  way  and  that  the  large  gap  in  scale  cannot  be  explained  in  terms  of 

3 1 

aliasing  effects  (I  would  accept  a  factor  of  S,  but  not  1000—  see  Gifford  )  or  other 
unexpected  turbulence  properties,  for  example,  an  anisotropic  eddy  shape  of  1000; 
l,  which  has  never  been  observed  in  the  laboratory  or  in  boundary  layer  turbulence. 
The  conclusion  then  seems  Inescapable  that  HICAT  spectra  are  predominantly  due 
to  gravity  waves.  For  this  reason,  all  attempts  to  explain  the  HICAT  spectra 


29.  Anderson,  A.  D,  (1957)  Free  air  turbulence,  J,  Meteorol.  1(477-494, 

30.  Pond,  S. ,  Steward,  H.W.,  and  Burling,  R.W.  (1963)  Turbulence  spectra  In 

the  wind  over  waves,  J,  Atm.  Scl,  20(319-321, 

31.  Gifford,  F,  (1959)  The  Interpretation  of  meteorological  spectra  and  correla¬ 

tions,  J.  Meteorol.  16(344-340. 
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on  the  basis  of  any  sort  of  turbulence  phenomena  must  be  regarded  as  Invalid  1 
A  quantitative  attempt  to  explain  the  HICAT  spectral  shape  In  terms  of  gravity  waves 
will  be  found  In  Dew  an, 


3.  THE  WORK  OF  HECK  AND  PANOFSKY 

The  starting  point  for  these  authors  Is  the  definition  of  eddy  dlffualvtty,  K^,  In 
terms  of  heat  flux 


% 


FLUX 

86 

7Z 


(14) 


that  Is, 


“h 


86/8Z 


(15) 


which  Is  Eq.  (5),  They  depart  from  the  direction  of  Lilly  et  al,  *5,  10  and 
23  -  . 

Zimmerman  et  al,  by  relating  Kjj  to  w'8'  directly  rather  than  to  t.  This  value 
of  w'61  Is  estimated  directly  from  the  HICAT  data  by  Integrating  the  cospeotrum  of 
w'  and  8',  that  is,  the  real  part  of  the  cross  power  spectal  density  of  6'  and  w'. 

In  order  to  remove  "random  effects  and  the  Influence  of  gravity  waves"  they  de¬ 
fined  0'w'  as  the  area  under  the  cospectrum  for  wavelengths  less  than  3000  m, 

Then  they  estimated  K„  from  Eq.  (15).  To  compare  these  values  of  with  those 
obtained  by  Lilly  et  al,  15,  15  they  computed  e  from  the  longitudinal  wind  component 
and  plotted  by  vs  log  e.  They  found  a  perfect  fit  with  =>  e/3Ng,  On  this 
basis  they  claimed  that  Kjj  was  of  order  104  cm2/s  and  with  F  —  0,01  In  Eq.  (13), 

Kg  -  0.  01  m2/s.  They  thus  fully  supported  the  results  given  by  Lilly  et  al.  l5,  10 
The  large  scale  cut  off  of  X  =  3000  m  exceeds  what  appears  to  be  the  outer 
length  by  a  factor  of  10,  The  choice  of  this  wavelength  Is,  therefore,  somewhat 
arbitrary.  Not  only  that,  but  It  conflicts  with  the  X  =  CIO  m  used  by  Lilly  et  al,,15’  10 
If  this  work  were  t.o  be  extended,  It  would  be  extremely  helpful  If,  In  addition  to  the 
cospectra,  the  quadrature  spectra  were  measured.  (This  is  the  Imaginary  part  of 
the  cross  spectrum. )  In  this  way  one  could  make  use  of  the  excellent  suggestion 


32. 


ewan,  E.M.  (1976)  Theoretical  Explanation  of  Spectral  Slopes  In  St 
spheric  Turbulence  Data  and  Implications  for  Vertical  Transport. 

TfffeTITTtTi'CTTTr'ATSTWBW.  - — 


33.  Welnstoek,  J,  (1976)  On  the  theory  of  turbulence  in  the  buoyancy  subrange  of 
stably  stratified  flows,  J,  of  the  Atms,  Set.  35SB34-049, 
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P  3  4 

which  was  made  by  Stewart  '  and  which  was  subsequently  demonstrated  by  Axford. 
This  suggestion  Involved  a  technique  for  differentiating  between  turbulence  and 
waves  In  the  atmosphere.  Defining  (3  to  be  the  phase  angle: 


<18) 


(where  Sw,0,  Is  the  cross  spectrum  between  potential  temperature  fluctuation  and 
vertical  velocity  fluctuation,  Im  designates  the  Imaginary  part  and  Re  the  real  part) 
they  found  that  one  could  distinguish  between  turbulence  and  waves  as  follows. 

When  the  flow  Is  mainly  turbulent  one  would  expect  0  to  be  In  the  range  4-45°  to 
■•45°  or  135°  to  226°.  For  flow  which  Is  essentially  wavellke,  /3  ■<  00°  ±  10°  (or 
270  ±  10°)s  and  coherence  values  greater  than  0,  8  are  usually  found  In  a  well- 
defined  wave  train,34  If  this  were  done,  a  less  arbitrary  criterion  than 
(X  <  3000  m)  might  be  obtained,  In  my  opinion  It  seems  quite  possible  that  the 
results  and  conclusions  of  Heck  and  Panofsky17  could  be  significantly  altered  by 
such  extensions  of  their  work, 


4.  THE  WORK  OF  ZIMMERMAN  AND  LOVING 


In  a  manner  similar  to  the  work  previously  cited,  these  authors  assumed  that 

the  HICAT  spectra  represented  Inertial  range  turbulence  and  that  one  could  use 

Eq.  (1)  for  the  purpose  of  estimating  e.  In  order  to  estimate  Kjj,  however,  they 

proceeded  quite  differently  and  avoided  the  contradiction  between  Eqs.  (1)  and  (3) 

that  has  been  discussed.  They  used  the  work  of  Heisenberg33  where  a  turbulent 

36 

dlffuslvlty  was  obtained  by  means  of  a  dimensional  argument,  namely 


W  ■  C1  rm 

ki 


dk 


(17)* 


There  la  an  element  of  arbitrariness  In  Eq,  (17)  tn  the  aense  that  any  K  of  the  form 


%  °  (constant) 


00 


4  Ik) 


s/2  k-f(s+iy2) 


d  k 


1/s 


or  combination  thereof  will  have  the  appropriate  dimensions,  Perhaps  the  best 
discussion  of  this  subject  will  be  found  In  Lin  et  al,  3 "»  38 

Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  31, 


IS 


Prom  Eq.  (1) 


KjjOtj)  -  a1/2  Cl  tl/3  kj4^3  .  (18) 

These  authors  then  addressed  the  question  of  what  value  to  use  for  k^.  For 

this  purpose  they  referenced  a  review  by  Phillips4  which  treated  the  "buoyancy 

39  40 

subrange"  of  turbulence.  In  particular  they  used  the  work  of  Lumley  and  Shur 
whose  theory  predicts  a  -3  dependence  of  the  spectrum  upon  k  (that  Is,  4  ~  k’3)  at 
scales  significantly  larger  than  the  Inertial  range  and  a  transition  wave  number, 
kB  given  by 

kB  *  c  N^2  .  (10) 

The  so-called  "buoyancy  length,"  was  given  by  *B  *“  and  they  estimated  It  to  be 
In  the  range  15  m  to  SI  m.  They  chose  C  ■  1  arbitrarily  In  order  to  calculate  this. 
They  then  chose  k^  a  kg.  It  should  be  noted,  however,  that  the  small  scale  resolu¬ 
tion  of  the  HICAT  data  was  close  to  30  m. 

Next,  In  order  to  use  Eq.  (18)  to  obtain  Kg  (k^)  from  «,  they  had  to  obtain 
values  for  a  and  Cj.  For  this  purpose  they  used  a  =  0.  5  (based  upon  published 

experimental  values).  For  C.  they  used  0.  B1  which  was  derived  In  an  appendix  by 
41  1  42 

Puo  and  Zimmerman.  This  appendix  used  data  published  by  Kellogg  and  used 

the  assumption  that  Eulerlen  and  Lagranglan  diffusion  were  approximately  equal. 

In  this  way  they  arrived  at  Kg  —  («/Ng)  (2.77)"1  which  Is,  for  present  purposes, 

Identical  to  Eq.  (13).* 

Three  specific  runs  from  the  HICAT  data  were  employed:  264  run  16,  233 

run  3,  and  280  run  10.  They  obtained  values  of  K„  ranging  from  2.  8  X  10®  cm2/s 

3  2  **  2  3 

to  37.  5X10  cm  /s  which  corresponded  to  a  range  of  e  from  24  cm  / s  to 

2  3  2 

92  cm  / a  .  The  largest  value  oft  was  listed  as  262  cm  /s.  These  should  be 

compared  to  values  found  by  Lilly  et  al,  where  e  (mean  dissipation  rates) 

93  2  3 

were  In  the  range  2.82  cm  ,/s  to  20,0  cm  /a  ;  and  Kg  values  given  In  tables  were 

jft 

ThlB  observation  was  originally  made  by  Dr,  It.  E.  Good  In  private  conversation. 

39,  Lumley,  J.L,  ( 1965)  Theoretical  Aspects  of  Research  on  Turbulence  In 

Stratified  Flows —Atmospheric  Turbulence  and  Radio-wave  Propagation, 
International  Colloquium,  Moscowsl05-112, 

40.  Shur,  G.  (1902)  Experimental  Studies  of  the  Energy  Spectrum  of  Atmospheric 

Turbulence,  Proceedings  of  Central  Aerologtcal  Observatory,  43. 

41.  Pao,  M.  and  Zimmerman,  S.  (1975)  Lagranglan  description  of  turbulent 

diffusion.  The  Natural  Stratosphere  of  1074.  CIAP  Monograph  1.  Final 
Report,  DoT-'T3T^'&T7(f- 1 23"- rt -124; - 

42,  Kellogg,  W.  (1956)  Diffusion  of  smoke  in  the  stratosphere,  Journal  of 

Meteorology  13:241. 
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in  the  range  of  2600  cm  /s  to  12,  200  cm  /a.  We  see  that  reasonable  men  can 

differ  by  two  orders  of  magnitude  when  It  comes  to  stratospheric  diffusion  estimates. 

When  estimates  as  obtained  by  Zimmerman  et  al  are  reduced  by  a  factor  of 

100,  they  fall  into  the  range  0.  3  to  3,  0  m2 /a  which  Is  the  same  range  reported  In 

the  literature  for  all  transport  mechanisms  for  stratospheric  transport. 

*3 

One  particularly  disturbing  aspect  of  their  treatment  is  that  the  k  dependence 
of  the  velocity  power  density  spectrum  (predicted  by  the  same  theory  that  they  used 
to  calculate  which  In  turn  was  used  in  Eq,  (IB)  to  obtain  K^j)  is  in  blatant  con¬ 
tradiction  to  the  HICAT  data  upon  which  they  base  their  estimates.  They  did  not 
overlook  this  fact  and  they  warned  the  reader  about  it,  They  did  not,  however, 
give  an  argument  as  to  why  this  would  not  invalidate  their  estimate  and,  therefore, 
we  are  left  with  this  as  an  unanawercd  question.  In  any  case,  if  their  results  were 
In  fact  of  the  correct  order  of  magnitude,  such  results  would  be  considered  aB 
evidence  that  turbulence  plays  a  significant  role  in  vertical  stratospheric  transport. 


S.  FURTHER  CRITICAL  REMARKS  APPLICABLE  TO  ALL  THE  ABOVE 
TREATMENTS 

S.1  111*  -5/3  Spectra  and  the  Determination  of« 

In  all  of  the  estimates  of  K„  described  to  this  point,  the  assumption  was  made 
that  the  HICAT  spectra  were  in  the  inertial  range.  The  assumption  seems,  how¬ 
ever,  to  be  not  valid  for  the  reasons  cited  previously.  On  the  other  hand,  the 
possibility  exists  that  in  spite  of  this  these  spectra  might  be  analogous  to  inertial 
range  spectra  to  some  degree  of  approximation. 

A  simple,  theoretical  explanation  of  the  HICAT  -5/3  spectra  was  proposed  by 

q  in  43 

Dewan  '  *  which  assumes  that  these  spectra  are  mostly  due  to  waves  and  that 

the  olope  Is  due  to  an  energy  cascade  caused  by  the  small  nonlinear  Interactions 
between  the  waves  of  various  scales,  This  wave  cascade  was  presumed  to  be  the 
Source  of  the  energy  which  eventually  finds  its  way  to  the  turbulence  cascades 
Inside  the  bltnl,  Thus,  In  fact,  If  this  theory  were  correct,  It  Implies  that  Indeed 
there  Is  a  physical  basis  for  the  existence  of  spectra  at  large  scales  which  are 


jfc  2 

S.  Zimmerman,  in  a  recent  private  communication,  revised  K«  down  to  0,1  m  /s 
on  the  basis  of  statistics  not  previously  considered.  The  conclusion  remains 
unaltered. 

43.  Dewan,  E.M.  (1979)  Mixing  in  Billow  Turbulence  and  Stratospheric  Eddy 
Diffusion.  Aroi,-flTT^WT7T17M7ri57r  K  * 


analogous  to  Inertial  range  spectra.  *  This  raises  the  next  question;  Would  the 
value  of  <  derived  from  such  spectra  be  even  approximately  equal  to  the  value  of 
e  found  In  a  valid  manner  (that  is,  from  an  actual  inertial  range  spectrum)? 

To  answer  this  question  we  consider  the  hypothetical  relation  between  the 
wave-cascade  spectrum,  S^(k),  and  the  inertial  turbulence -cascade  spectrum, 
ST(k)t 

5^y(k)  *  a w*w^^ 

from  Dew  an0,  10'  43  and 

ST(k)  «  aT  «T2/3  kT‘6/3  (21) 

where  subscripts  W  and  T  are  used  to  designate  wave  and  turbulence  quantities. 
While  aT  is  known,  all  that  one  knows  about  a^  Is  that  it  is  of  order  unity, 

According  to  Bond,  we  do  know  a  little  more.  He  haB  demonstrated  that  con¬ 
stants  like  aw  are  either  greater  than  S  or  less  than  1/S  in  1/8  of  all  cases  studied 
by  him  and  that  the  probability  Is  1/10  that  it  be  greater  than  10  or  less  than  10, 
etc,  For  convenience  In  what  follows,  we  shall  set  aw  *  aT  and  try  to  relate  eT 

Let  us  assume  that  the  flux  of  energy  from  large  scale  (over  10  km)  to  small 
scale  (less  than  1  m)  Is  conserved.  In  other  wordB,  we  assume  that  there  are  no 
energy  sources  or  sinks  at  Intermediate  scale,  This  would  be  a  dubious  assumption 
in  the  case  of  boundary  layer  turbulence  since  heating  of  the  ground  can  cause 
buoyancy  driven  turbulence  of  Intermediate  scales.  We  assume,  however,  that 
Eqs.  (20)  and  (21)  represent  parts  of  a  conservative  cascade.  This  leads  directly 
to  the  following  expression  of  flux  conservation  (down  the  scale); 


One  other  possibly  viable  suggestion  has  been  proposed  to  explain  the  type  of 
spectra  seen  In  the  HICAT  results.  ">  44  ThlB  explanation  involves  the  ''two 
dimensional"  turbulence  cascade  which  according  to  Kralchnan,  45  gives  a  k'5/3 
spectrum  (one-dimensional)  associated  with  a  cascade  going  In  the  reverse  direc¬ 
tion)  that  is,  from  small  to  large  scales.  The  explanation  leaves  unanswered 
three  questions;  Why  does  the  two-dimensional  turbulence  exist  In  layers  of 
order  1  km  In  thickness?22  (This  Is  easily  explained  in  terms  of  trapped  gravity 
waves  on  the  other  hand.)  Also,  what  is  the  source  of  energy  at  small  scales 
(leas  than  50  m)  ?  If  the  answer  is  "turbulence,"  then,  what  is  the  source  for  the 
turbulence  energy? 

44.  Gage,  K,  S.  (1970)  Evidence  for  a  k’6/3  low  Inertial  range  In  mcsoscale  two- 

dimensional  turbulence,  J,  Atm,  Scl,  30 1 1050-1054. 

45.  Kralchnan,  R.H.  (1067)  Physics  of  Fluids  10;  1417. 

40.  Bond,  W.  N.  (1929)  Phil.  Mag.  7;722. 
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where  B,^  la  the  energy  dissipation  rale  in  the  wave  oaaoada  par  unit  volume, 
the  aama  quantity  for  the  turbulence  caacade,  and  VW  and  VT  the  volume  ocoupied 
by  the  wavea  and  turbulence  reapecttvely.  That  theae  volumea  oan  differ  slgnlfl- 
cantly  waa  demonatrated  by  the  obaervatlona  of  Wooda,  *  and  Woods  and  Wiley.  ® 
They  saw  Internal  wavea  in  the  upper  ocean  perturbing  relatively  rare  and  thin 
turbulent  layera  and  yielding  their  energy  (originally  apread  over  both  laminar  and 
turbulent  layera  alike)  to  the  turbulence  caacadea  via  the  Kelvln-Helmholta  Insta¬ 
bility,  Now 

m  p«w  and  DT  ■  p*T  (23) 

where  p  la  the  maaa  density.  Thus,  from  Eq.  (22 ) 

<P«W>VW  »  <p«T)VT  (24) 


Thus,  If  Vw  »  VT,  which  la  to  be  expected,  then  eT  »  and  KH  would  have  to 
be  altered  accordingly  (assuming  that  -  c  were  valid). 

An  Important  but  unanswered  question  In  connection  with  the  above  argument 
Is  "are  the  HICAT  spectra  due  mostly  to  waves  ?"  This  could  be  determined  by 
means  of  the  appropriate  cross  spectra,  but  such  has  not  yet  been  done  In  the 
literature,  Due  account  would  also  have  to  be  taken  of  aliasing  effects  In  such  a 
study,  But  In  spite  of  such  unanswered  questions,  It  Is  now  posaible  to  state  with 
certainty  that  It  la  Indeed  possible  to  make  a  significant  error  In  the  estimation  of 
e  directly  from  -5/3  spectra  when  In  fact  the  latter  represent  waves  Instead  of 
turbulence  as  Is  moat  likely  the  case  for  HICAT  data.  This  would  Imply,  for 

2  I  a  1 

example,  that,  assuming  Kj}  -  «/NB  were  valid,  the  resultB  of  Lilly  et  al,  ' 
are  too  low  by  the  factor  (V^/V^). 


&  n  in 

These  considerations  were  not  Included  In  Dewan.  •  1  1 


S.2  Extreme  Inhomogeneitv  end  Ite  Effect  lij>on  the  Keiatioiuhip 
Between  K|(  and  Kj| 

In  the  following  It  will  be  shown  that  there*  are  conditions  such  that  Kg  would 
have  no  influence  upon  the  value  of  Kg  and  that  such  conditions  are  not  at  all  un¬ 
likely  In  the  case  of  stratospheric  turbulence  (or  stratified  turbulence  in  general). 
As  in  Section  S,  1,  all  the  above  methods  to  estimate  Kg  will  be  affected  by  the 
criticisms  raised  there,  Eq.  (13)  la  repeated  here  for  convenience, 

Kg-EKg  .  (13) 

In  other  words,  one  multiplies  the  eddy  diffusion  coefficient  corresponding  to  tho 
turbulence  effects  within  an  actively  mixing  layer  by  the  fraction  of  the  vertical 
dimension  which  is  turbulent  in  order  to  estimate  the  bulk  transport  In  that  direc¬ 
tion  in  terms  of  a  diffusion  parameter.  One  "dilutes"  Kg  to  estimate 
almost  anyone  this  sounds  reasonable  at  first.  Eq.  (13),  however,  is  not  always 
valid.  The  layered  turbulence  structure  (blinl)  In  the  stratosphere  has  already 
been  described;  and  these  rare,  active,  mixing  layers  separated  by  laminar  flow 
bear  no  resemblenco  whatever  to  homogeneous  turbulence.  Usually  homogeneity 
assumptions  are  made  In  turbulence  theory,  but  it  should  be  clear  that  in  the 
present  case  any  such  assumption  would  be  manifestly  invalid.  We  shall  examine 
below  some  special  examples  where  Eq.  (13)  la  clearly  Invalid  so  that  the  main 
difficulties  associated  with  It  will  be  put  in  evidence. 

First,  It  Is  necessary  to  define  Kg  as  explicitly  as  possible.  Figure  1  shows 
a  slab  of  atmosphere  much  thicker  than  a  turbulent  layer.  Let  a  scalar  constituent 
(temperature,  neutrally  buoyant  trace  gss,  etc. )  be  given  by  *(2)  and  let  A  and  B 
be  points  at  the  top  and  bottom  of  the  slab  as  Indicated.  We  define  Kg  by 

K  m  FLUX  (from  A  to  B) 

B  <8?/9Z) 

where,  from  A  to  B  the  gradient  9V/QZ  can  be  taken  as  constant, 

We  now  consider  the  case  whore  there  is  a  single  active  layer  of  turbulence 
located  at  altitude  C,  and  this  layer  is  presumed  to  be  permanently  fixed  for  all 
time,  For  example,  let  the  distance  from  A  to  B  be  10  km  and  the  layer  thickness 
be  200  m.  According  to  Eq.  (13),  Kg  -  Kg  (200  m/10*  m).  Is  this  correct? 
Consider  Eq.  (26).  Since  there  is  obviously  no  flux  from  point  A  to  point  B  taking 
place  (Figure  1),  one  must  conclude  that  in  the  case  considered,  Kg  =  0.  (Molec¬ 
ular  transport  is,  of  course,  being  omitted  from  the  discussion  there.) 

*  12 
The  first  part  of  Section  Pi.  2  is  a  brief  summary  of  part  of  Dewan. 
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Figure  1.  Single  Turbulent  Layer  In  Slab  of  Fluid 
Extending  Between  Altitudes  A  and  B,  The  & 
denotefl  the  mixing  ratio  of  pollutant 


The  above  situation  resembles  the  ease  In  electrical  conductivity  where  one 
has  a  "sandwich"  consisting  of  a  horizontal  slab  of  Insulation  above  and  below  a 
horizontal  slab  of  conductor  as  shown  In  Figure  2,  In  this  configuration  there  la 
no  way  that  current  can  flow  vertically  as  Indicated. 


ZERO  CURRENT 
THROUGH  "SANDWICH” 


Figure  2.  Insulator-Conductor-Insulator  Sandwich 
Analogy  for  Transport  by  Layered  Turbulence 


Next,  consider  Figure  3  which  depicts  a  more  realistic  situation  where  the 
single  fixed  layer  is  replaced  by  an  ensemble  of  randomly  spaced  layers.  Their 
thicknesses  are  random  as  Indicated  as  well  as  their  spacing,  Total  mixing  is 
presumed  to  take  place  within  the  turbulent  layers,  It  we  assume,  as  before,  that 
these  layers  remain  at  fixed  altitudes  for  the  entire  duration  of  the  experiment, 
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then,  again,  we  must  conclude  that  Kg  «  0  (since  again  there  Is  no  transport  from 
point  A  to  point  B).  On  the  other  hand,  let  us  next  suppose  that  the  layers  do  not 
maintain  a  constant  configuration  but,  Instead,  have  a  finite  duration  and  are 
replaced  by  other  layers  at  random  heights  and  thicknesses.  After  a  sufficient 
length  of  time  has  elapsed,  all  altitudes  will  have  been  covered  by  turbulent  layers 
several  times.  Thus,  when  enough  time  has  passed  there  would  Indeed  be  a  flux 
from  A  to  B.  This  type  of  flux  has  been  studied  In  detail  In  Dewan11  by  means  of 
computer  simulations,  The  flow  over  short  times  will  be  very  Irregulari  but  over 
long  times  and  with  stationary  statistical  behavior  of  layer  formation,  the  average 
flow  becomes  constant,  Thus  Kg  could  be  directly  "measured"  by  means  of 
Eq.  (20).  12 

Let.  At  be  the  time  Interval  between  different  layer  configurations.  In  other 

words,  the  time  between  the  commencement  of  one  entire  layer  ensemhle  (Figure  3) 

and  the  event  when  It  In  replaced  by  a  new  one  Is  At,  If  At  were  doubled,  what 

would  happen  to  the  flux  going  from  A  to  B?  It  would,  of  course,  be  halved,  and 

hence  KB  would  be  reduced  by  a  factor  of  2.  From  this  we  see  that  Kg  ~  (At)"1. 

In  what  way  will  Kg  depend  upon  A,  the  average  layer  thickness?  Also,  how  would 

Kt,  depend  upon  F,  the  fraction  of  vertical  dimension  turbulent?  Since  these 
a  1112 
questions  have  been  treated  elsewhere,  and  since  the  answer  can  be  obtained 

by  using  Eq.  (28)  directly,  we  merely  state  the  final  result,  10> 


K 


B 


F  A 
”21 


(27) 
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Here  >  «*  (1/10).  The  most  Important  thing  to  notice  In  Eq.  (27*)  1b  that  Kg  (or  e) 
is  not  relevant.  Only  the  layer  thickness,  A,  "cycle  time,"  At,  and  F  play  any 
role  In  the  flux  In  Figure  3  (A  to  B).  The  assumption  that  nearly  total  mixing  takes 
place  la  the  reason  that  Kg  drops  out  or  the  argument  in  this  inhomogeneous  situa¬ 
tion.  Thus,  under  such  conditions  one  must  rule  out  Eq,  (13)  for  Kg. 

The  parameter,  Kg,  given  by  Eq,  (27)  is  not  an  eddy  diffusion  coefficient  of 
the  usual  type,  Eddy  diffusion  la  usually  associated  with  the  definition, 


K-v'l'  (28) 

where  v1  and  1*  are  the  velocity  and  length  scales  of  homogeneous  turbulence.  To 
avoid  confusion,  therefore,  perhaps  one  should  call  Kg  In  Eq.  (27)  the  "stratified 
turbulenoe  diffusion  parameter." 

Does  nearly  total  mixing  take  place  in  the  stratosphere  ?  This,  of  course,  is 

one  of  the  most  Important  questions  in  connection  with  the  use  of  Eq.  (27).  Mantis 
47 

and  Peppln  measured  temperature  profiles  in  the  upper  troposphere  and  lower 
stratosphere  by  means  of  balloon-borne  sensors.  Their  observations  are  consis¬ 
tent  with  the  hypothesis  that  nearly  total  mixing  occurs,  because  they  found  nearly 
adiabatic  lapse  rates  over  regions  of  order  200  m  thick  which  numerically  Is  the 
same  thickness  to  be  expected  of  typical  turbulence  layers  (a  "coincidence"  not  to 
be  taken  lightly).  The  only  alternative  explanation  I  can  find  Is  that  the  adiabatic 
lapse  rate  regions  may  actually  be  unrelated  to  turbulence  but  be  actually  due  to 
Intrusions  of  layers  of  air  (at  particular  altitudes)  which  are  premixed.  *  Such  an 
hypothesis  while  it  would  explain  the  observations  would  leave  unanswered  the 
questions  of  how  and  where  such  intrusion  layerB  could  be  formed  In  the  strato¬ 
sphere. 

2  ^ 

Observations  In  the  upper  ocean  ’  show  temperature  profiles  which  are 
shaped  like  "steps"  (that  Is,  have  regions  where  temperature  Is  constant  with 

respect  to  depth).  The  else  of  these  steps  (in  complete  analogy  with  the  observe- 

47 

tlons  of  Mantis  and  Peppln  )  Is  the  same  as  the  typical  size  of  n  turbulent  layer 
in  the  ocean.  Since  there  Is  some  analogy  between  the  dynamics  of  the  upper  ocean 
and  stratosphere,  one  can  regard  this  observation  as  Indirect  evidence  of  the  total 
mixing  assumption.  On  the  other  hand,  the  "Intrusion  layer"  phenomenon  la  known 
to  occur  In  the  ocean. 

While  there  exists  published  evidence  which  Is  consistent  with  the  nearly  total 
mixing  assumption,  and  while  there  Is  numerical  agreement  between  the  vertical 


% 

Koop  In  private  communication. 

47.  Mantis,  M.  T.  and  Pepin,  T.J,  (1971)  Vertical  temperature  structure  of  the 
free  atmosphere  at  mesoscale,  J.  Geophys,  Res,  20i 8 02 1  -8028. 
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thickness  of  mixed  layers  and  the  thicknesses  of  observed  turbulent  layers,  this 

issue  is  not  yet  settled.  In  view  of  this,  it  is  Important  to  consider  reasonable 

experiments  which  can  be  performed  to  obtain  more  definitive  evidence.  An 

especially  attractive  approach  seems  to  be  the  technique  used  by  Browning  and 
46  49 

Watkins  and  Browning.  Under  certain  conditions  the  troposphere  exhibits  the 
same  sort  of  bltnt  as  those  seen  In  the  stratosphere.  These  are  called  CAT  (clear 
air  turbulence)  and  appear  on  radar  in  a  very  well-defined  manner.  These  authors 
used  radar  observations  on  tropospheric  CAT  layers  in  conjunction  with  balloon 
observations  of  temperature  and  wind  shear.  They  were  able  to  measure  Richardson 
number  profiles  as  a  function  of  time  by  means  of  repeated  balloon  launchings.  In 
Figure  3  of  Browning  and  Watkins  the  vertical  gradient  of  potential  temperature, 
the  vortical  shear  of  the  horizontal  winds,  and  the  vertical  profile  of  R^  are  plotted 
before  and  after  the  occurrence  of  a  billow  event,  This  figure  clearly  shows  that 
some  mixing  occurs,  but  if  this  were  typical,  it  would  indicate  that  the  mixing  is 
significantly  far  from  total.  Unfortunately,  the  resolution  of  the  balloon  measure¬ 
ments  Is  about  200  m,  and  higher  resolution  would  be  very  helpful.  Also,  if  billows 
typically  did  not  cause  a  large  amount  of  mixing,  there  remains  the  question  of 
whether  or  not  a  series  of  billows  such  an  those  seen  in  references  2  and  3  could 
in  time  cause  nearly  adiabatic  lapse  rates.  If  this  were  tD  be  the  case,  then  one 
would  have  to  estimate  At  with  this  in  mind,  and  it  appears  likely  that  At  would  then 
be  significantly  more  difficult  to  estimate  than  otherwise  (set*  below).  The  effect 
would  be  to  increase  At  and  decrease  Kg. 

Alternatively,  If  much  less  than  total  mixing  were  to  occur,  Eq,  (27)  could  be 
modified  to  take  this  Into  account,  Such  possible  modifications  will  be  discussed 
elsewhere!  however,  for  sufficiently  small  mixing,  it  has  been  shown0,  *3  that 
Eq.  (27)  for  Kg  would  have  to  be  replaced  by  Eq.  (13).  Thus,  there  are  conditions 

where  Eq.  (13)  is  actually  valid.  To  use  Eq.  (13),  however,  as  Lilly  et  al10’ 

1 7 

and  Heck  et  al  have  done,  leaves  us  with  unanswered  questions  pertaining  to  its 
validity.  First,  one  must  obtain  valid  experimental  evidence  of  the  actual  degree 
of  mixing  In  stratospheric  turbulent  layerB,  In  u  word,  perhaps  one  of  the  moat 

Important  single  experiments  to  perform  at  this  time  regarding  stratified  turbulence 

47 

would  be  a  repetition  of  the  work  of  Browning  and  Watkins  with  higher  resolution 
and  over  longer  periods  of  continuous  observations, 


48.  Browning,  K.A,  and  Watkins,  C'.D,  (1D70)  Observations  of  clear  air  turbulence 

by  high  power  radar,  Nature  227i260-2fl3, 

49,  Browning,  K,  A,  (1971)  Structure  of  the  atmosphere  In  the  vicinity  of  large 

amplitude  Kelvtn-Helmholtz  billows,  Roy,  Met.  Soc,  Quart,  J.  97i283-298. 
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6.  THE  WORK  OF  ROSENBERG  AND  DEWAN 


O 

Equation  (27)  was  first  derived  and  applied  by  Rosenberg  and  Dewan.  It  was 
applied  to  a  data  base  consisting  of  200  vertical  profiles  of  horizontal  winds  orig¬ 
inally  obtained  by  Miller  et  al80  of  NASA,  These  wind  profiles  were  obtained  by 
sending  up  rockets  which  left  behind  long  trails  of  Bmoke  which  were  approximately 
vertical.  These  smoke  trails  were  then  photographed  from  three  widely  spaced 
locations  In  simultaneous  time-lapsed  fashion,  The  series  of  photographs  were 
used  to  estimate  wind  by  means  of  trlangulatlom  The  location  of  points  along  a 
trail  could  be  obtained  In  three  dimensions  as  a  function  of  time,  and  this  gave 
the  horizontal  wind  profile.  The  altitude  range  studied  was  12  to  18  km,  and,  thus, 
It  Included  both  upper  troposphere  and  lower  stratosphere.  Miller  and  Henry 

KQ 

et  al,  provided  the  data  at  a  vertical  resolution  of  25  m, 

Q 

In  order  to  estimate  the  turbulent  transport  effects,  Rosenberg  and  Dewan 
first  obtained  vertical  profiles  of  the  Richardson  number,  The  main  Idea  was  to 
ascertain  which  altitudes  had  R^  <  1/4,  and,  hence,  estimate  the  probability  of 
occurrence  and  probable  thicknesses  (P  and  A  In  Eq,  (27))  of  the  turbulent  layers, 
Unfortunately,  no  simultaneous  temperature  profiles  were  available  In  connection 
with  the  wind  profiles,  and,  hence,  in  order  to  obtain  the  Richardson  numbers 
(which,  of  course,  depend  on  the  potential  temperature  gradients),  we  had  to  re¬ 
sort  to  a  model  atmosphere.  This  raises  the  as  yet  unanswered  question!  How 
much  would  the  results  be  altered  If  actual  rather  than  assumed  temperature  pro¬ 
files  were  employed?  Our  laboratory  has  conducted  experiments  which  should 
help  to  provide  an  answer  to  this  question  In  the  future. 

A  second  unanswered  question  relating  to  this  work  tnvolves  the  fact  that 
sometimes  R^  could  go  below  1/4  and  then  return  subsequently  to  a  value  above  1/4 

without  an  intervening  billow  event,  This  effect  would  have  to  be  Included  In  accu- 

49 

rate  assessments  of  Kg.  The  work  of  Browning  shows  that  this  does  Indeed 
happen  sometimes. 

This  work  took  Into  account  the  effect  of  turbulent  spreading.  As  has  already 

2 

been  mentioned,  the  conclusion  was  that  Kg  —  0,3  m  / s  which  would  make  turbu¬ 
lence  of  the  small  scale  variety  one  of  the  significant  agents  In  the  vertical  trans¬ 
port  proceas  within  the  stratosphere.  Total  mixing  within  the  layers  was  assumed, 
of  course,  and  this,  as  was  already  emphasized,  must  be  further  tested.  But 

g 

there  la  yet  one  more  assumption  made,  not  only  by  Rosenberg  and  Dewan,  but 
by  all  the  previously  cited  authors,  which  tf  not  correct,  could  wreak  havoc  with 


30.  Miller,  R.W.,  Henry,  R.M.,  and  Rowe,  M.G,  (1935)  Wind  Velocity  Profiles 
Measured  by  the  Smoke-Trail  Method  at  Wallops  Island,  Virginia  NASA 

Tw&mrtimttmr;  see  sIboTeasa  rsrp-w'araBtfBy:  — 
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all  Kg  estimates.  This  Is  the  assumption  that  turbulent  layers  form  at  random 
altitudes. 

The  postulate  of  the  randomness  of  altitude  for  turbulent  layer  formation  has 
the  following  justifications!  (1)  It  Is  not  in  contradiction  with  known  data,  and 
(2)  It  Is  the  simplest  hypothesis  to  make  in  the  face  of  our  present  Ignorance.  On 
the  other  hand,  It  still  remains  entirely  possible  that  turbulent  layers  have 
"preferred  altitudes"  over  long  durations  of  time  relative  to  the  time  history  of 
layer  creation  and  subsequent  decay.  Such  a  possibility  Is,  In  fact,  Implicit  in  the 
previously  cited  theory,  Dewan,  9’  10,  *3  for  the  explanation  of  the  HICAT  data.  In 
the  latter  theory  It  Is  assumed  that  project  HICAT  measured  trapped  gravity  waves. 
Trapped  gravity  waves  could  certainly  give  rise  to  preferred  altitude  regions  of 
turbulent  layer  formations. 

The  best  technique  for  testing  the  random  altitude  assumption  would  be  one 
which  uses  high  powered  radar  to  detect  stratospheric  turbulence  at  one  geographic 
location  over  extended  periods  of  time.  Such  observations  have  been  made  by 
Crane,  28,  61  VanZandt  et  al,  52  Woodman,  83  and  Watkins.  Unfortunately,  in 
all  but  one  of  these  observations  (Woodman's,  with  resolution  of  150  m)  the  highest 
resolution  Is  of  order  1  km.  Since  the  expected  layer  thickness  Is  of  order  100  m, 
the  resolution  must  be,  In  general,  greatly  Improved  and  more  extensive  observa¬ 
tions  must  be  made.  Slnoe  this  assumption  of  randomness  Is  crucial  for  the  theory 
(that  Is,  Kg  could  °  0  as  we  have  already  seen  In  Figure  1),  such  experiments  have 
very  high  priority  If  not  the  highest  priority  with  regard  to  turbulence  transport. 

In  Rosenberg  and  Dewan8  the  symbol  "L"  was  used  for  "layer  thickness." 

Later  on  It  became  evident  that,  In  actual  fact,  L  was  a  "half-thickness, 11  that  la, 

L  #  (A/2).  In  view  of  this,  we  found  that  A  typically  was  of  order  200  m.  Turbu- 
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lent  layers  of  such  size  were  measured  by  Barat,  Crane,  ’  and  Cadet,  by 
means  of  In-sltu  measurements,  These  measurements  lend  credibility  to  the  value 
of  A  which  were  reported  In  Rosenberg  and  Dewan,  even  though  we  subsequently 
learned  that,  at  20  m  resolution  there  was  a  large  error  In  velocity  (of  order 
0.1  m/s  to  1.0  m/s).  Further  assurance  came  from  the  fact  that  an  estimate  of 
Tv  based  on  a  smooth  velocity  profile  of  100  m  resolution  resulted  In 

22 

I  Ignore  here  the  observation  by  Crane  that  there  seems  to  be  a  "persistent  layer 
at  the  tropopause, "  (compare  his  Figure  32).  The  reason  is  that  we  are  consider¬ 
ing  transport  throughout  the  body  of  the  stratosphere  and  not  the  boundary  effects. 

51.  Crane,  R.K.  ( 1  ©80)  A  review  of  radar  observations  of  turbulence  In  the  lower 

stratosphere,  Radio  Science  1  St  1 77 -1 93 . 
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profiles  of  refractlvlty  turbulence  structure  constant,  Radio  Set.  13t8 19-82©. 
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54.  Watkins,  B.  (1981)  Radio  Science  (In  press). 
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2  2 

Kg  »  0,  24  m  /a  Instead  of  0,  34  m  /s  for  26  m  resolution.  Such  a  small  reduction 
In  Kg  did  not  alter  our  conclusion,  and  the  100  m  spacing  greatly  diminished  the 
spurious  effects  due  to  the  "error"  In  velocity  previously  mentioned.  Further 
treatment  of  this  situation  will  be  given  elsewhere. 

In  Rosenberg  and  Dewan,  (Figure  5  In  that  paper)  we  found  that,  at  25  m 
resolution  and  100  m  resolution,  the  fraction  F  of  turbulence  Is  of  the  order  of  a 
few  percent  (for  L  2  25  m),  This  fraction  of  turbulence  Is  numerically  of  the  same 
general  size  found  from  the  HICAT  data,  It  Is  tempting  to  regard  this  as  a  form  of 
mutual  confirmation}  however,  caution  Is  needed  here.  The  H1CAT  data  seem  to 
have  been  primarily  associated  with  gravity  waves.  In  addition,  they  were  taken 
over  nearly  horizontal  trajectories.  In  contrast,  the  F  In  our  work  Involved 
essentially  vertical  trajectories  and  involved  "potential  turbulence,"  that  Is, 
regions  where  <  1/4  together  with  the  estimated  effects  of  spreading,  That  F 
Is  about  equal  In  the  two  cases  may  thus  be  due  to  coincidence.  In  any  case,  the 
F  derived  by  means  of  vertical  profiles  Is  the  only  one  that  is  relevant  so  far  as 
vertical  Kg  Is  concerned  (Eq,  (27)). 

Another  parameter  In  Eq.  (27)  la  At.  In  Rosenberg  and  Dewan8  the  value  of 
At  was  estimated  from  the  previously  cited  work  of  Browning.  We  found  that, 
on  average,  3000  sec  would  elapse  between  the  time  R(  <  1/4,  according  to  the 
balloon  observations  and  the  billow  event  (turbulence)  according  to  the  radar  ob¬ 
servation.  Since  these  observations  were  carried  out  In  the  troposphere,  It  was, 
of  course,  necessary  to  make  allowance  for  the  change  of  the  dynamical  situation 
In  the  stratosphere,  and,  hence,  At  there  was  taken  as  1500  sec  (see  the  report 
for  details  of  the  argument),  I  feel  that  this  At  Is  accurate  enough  to  estlmato  Kg 
"within  one-half  an  order  of  magnitude"  1/  the  assumptions  are  correct.  However, 
only  after  the  key  experiments  where  degree  of  mixing  and  randomness  of  altitude 
for  turbulence  have  been  performed  would  It  make  any  sense  at  all  to  concentrate 
on  the  task  of  obtaining  more  accurate  estimates  for  At. 

One  more  critical  remark  should  be  made  regarding  the  estimate  of  Kg  In 
Rosenberg  and  Dewan.  8  It  was  based  entirely  upon  data  obtained  from  a  single 
geographic  location,  namely,  Wallops  Island,  Virginia,  In  addition,  the  data  were 
obtained  only  on  clear  days,  It  would  be  much  better  If  at  least  the  geographical 
effects  could  be  Included  since,  as  Lilly  et  al,  15,  have  demonstrated,  velocity 
fluctuations  change  significantly  with  changes  In  the  shape  of  the  terrain  beneath. 


7.  CONCLUSION 

We  have  seen  that  Kn,  the  bulk  vertical  transport  parameter  for  small  scale 
U  1  1  B  fi 

turbulence  has  been  estimated  to  bo  of  order  0.  01,  *  '  0. 10,  and  1. 0  or  even 
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The  question  Is  which,  1/  any,  of  these  estimates  Is  the  one  to  be  taken 

.23 


higher,  2S* 

seriously?  The  higher  values  of  Zimmerman  and  Loving”*  and  Rosenberg  and 

g 

Dewan  would  Imply  that  turbulence  plans  a  significant  role  in  global  vertical  trans* 

IK  1  fl 

port.  The  value  published  by  Lilly  et  al,  ’  would  Imply  an  Insignificant  role. 

We  have  seen  that,  without  exception,  all  these  estimates  Involve  assumptions 
and  unanswered  questions  which  can  only  be  tested  and  resolved  by  means  of 
experiments  that  have  yet  to  be  adequately  performed.  In  my  opinion  It  Is  clear 
that,  on  the  basis  of  the  wide  range  of  these  Independent  estimates,  and  the 
unanswered  questions  associated  with  their  derivations  the  most  valid  conclusion 
Is  that  at  present  no  available  value  of  Kg  Is  to  be  regarded  as  being  conclusive 
or  reliable,  All  we  really  know  Is  that  it  cannot  be  larger  than  the  values  obtained 
from  the  fallout  of  tracers  from  the  stratosphere  which  represent  all  the  processes 
operating  In  concert, 

But  now  let  us  bring  In  one  more  element,  Mahlman  and  Moxlm  In  a  paper 
entitled,  "Tracer  Simulation  Using  a  Global  General  Circulation  Modeli  Results 
From  a  Mid-latitude  Instantaneous  Source  Experiment"  discussed  the  significance 
of  the  role  of  what  they  called  "vertical  subgrld-scale  diffusion,"  (p.  1349), 18  In 
order  to  obtain  an  upper  bound  on  Kg  they  Inserted  values  of  0. 1  m2/s  and  0.  8  m2/s 
Into  their  "Global  Circulation  Model,"  They  state  their  results  as  follows!  "In 
both  cases,  the  tracer  transport  from  the  stratosphere  to  the  troposphere  was 
drastically  overestimated  compared  with  observed  behavior  of  radioactive 
tracers  , , ,  ,  The  above  result  suggests  that  subgrld-scale  motions  are  consider¬ 
ably  less  Important  than  large-scale  motions  In  affecting  stratospheric  vertical 
tracer  transfer,  This  Inference  is  strengthened  by  analysis  of  Project  HICAT 
spectra  (for  example,  Lilly  et  al, )."  18 

As  we  have  seen,  there  are  too  many  unanswered  questions  In  the  work  cited  by 
18 

Mahlman  and  Moxlm  to  lend  much  strength  to  their  Inference!  and  perhaps  It  Is 
better  to  regard  their  result  as  standing  oolely  on  Its  own  merits.  In  any  case, 
they  rightly  consider  It  to  be  desirable  to  have  an  Independent  way  to  ascertain  the 
value  of  Kg.  This  would,  therefore,  reinforce  the  need  to  perform  the  key  experi¬ 
ments  mentioned  In  the  text  regarding  degree  of  mixing  and  randomness  of  altitude 
of  formation.  If  Indeed  the<upper  limit  of  Kfi  Inferred  by  Mahlman  and  Moxlm,  18 
eventually  received  valid  Independent  support,  then  the  role  of  small  scale  turbu¬ 
lence  would  be  at  last  established.  On  the  other  hand,  It  would  still  be  Important 
to  obtain  more  than  the  upper  bound.  The  value  of  Kg  remains  Important  In  the 
context  of  the  chemistry  of  the  stratosphere. 


But,  us  was  mentioned,  this  has  recently  been  revised  downward  to  0, 1  m2/s  In  a 
private  communication. 
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Finally,  I  would  like  to  state  that  our  Ignorance  about  Kg  Is  much  larger  than 
may  have  seemed  possible  and  that  more  work  needs  to  be  done  In  view  of  the  Im¬ 
portance  of  the  problem. 
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